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ABSTRACT 

At least one massive binary system containing an energetic pulsar, PSR B1259- 
63/SS2883, has been recently detected in the TeV 7-rays by the HESS telescopes. 
These 7-rays are likely produced by particles accelerated in the vicinity of the pulsar 
and/or at the pulsar wind shock, in comptonization of soft radiation from the massive 
star. However, the process of 7-ray production in such systems can be quite compli- 
cated due to the anisotropy of the radiation field, complex structure of the pulsar 
wind termination shock and possible absorption of produced 7-rays which might initi- 
ate leptonic cascades. In this paper we consider in detail all these effects. We calculate 
the 7-ray light curves and spectra for different geometries of the binary system PSR 
B1259-63/SS2883 and compare them with the TeV 7-ray observations. We conclude 
that the leptonic IC model, which takes into account the complex structure of the 
pulsar wind shock due to the aspherical wind of the massive star, can explain the 
details of the observed 7-ray light curve. 

Key words: binary systems - pulsars: PSR B1259-63 - radiation mechanisms: non- 
thermal - gamma-rays: theory 



1 INTRODUCTION 

PSR B1259-63/SS 2883 is the best known massive binary 
system containing a classical radio pulsar in a highly ec- 
centric orbit around the Be type star. The pulsar creates 
strong relativistic wind which prevents accretion from the 
wind of the massive star. So then, it is expected that high 
energy processes, which are characteristic for the isolated 
radio (7-ray) pulsars, can also occur in the inner pulsar 
magnetosphere of PSR B1259-63. Due to the short rota- 
tional period, equal to 47 ms, and strong surface magnetic 
field, the pulsar produces a strong relativistic wind which, 
in contrary to isolated pulsars, interacts with the wind from 
the massive star. This interaction process is especially im- 
portant at the periastron passage when the pulsar is im- 
mersed in the inhomogeneous massive star wind. As a result 
of this interaction process, a relativistic shock in the pulsar 
wind appears. Particles can be accelerated on the shock to 
~ TeV energies. Therefore, it has been suspected, that the 
binary system PSR 1259-63/SS 2883 might produce observ- 
able high energy 7-ray radiation. This expectation has been 
confirmed during the last periastron passage of the pulsar 
in 2004, when the HESS detectors discovered the TeV 7-ray 
emission from this binary system (Aharonian et al. 2005). 

The production of X-ray and 7-ray radiation in PSR 
B 1259-63 system in the model of interaction of the pulsar 



and stellar winds was considered for the first time by Grove 
et al. (1997), and Tavani et al. (1997). In those works e 
pairs, escaping from the pulsar magnetosphere, are acceler- 
ated additionally at the pulsar wind termination shock and 
interact with the magnetic field and thermal radiation of the 
massive star. Moreover, the conditions for propagation and 
cooling of e + e~ pairs have been considered as a function 
of the orbital phase. The produced X-ray light curves were 
in good agreement with the observed X-ray flux variability. 
The photon spectra were fitted under the assumption that 
the massive star equatorial wind is inclined to the orbital 
plain with the inclination angle i > 25°. According to these 
models, the Lorentz factors of e pairs in the pulsar wind 
should be of the order of j e ~ 10 6 , while the wind mag- 
netization parameter, which determines the magnetic field 
strength at the shock, should be a ~ 10 _1 — 10 -2 (Tavani 
et al. 1996). 

A similar model for this binary system has been also 
considered by Kirk, Ball & Skjaeraasen (1999). These au- 
thors argued for the first time that also the TeV 7-rays 
should be produced in PSR 1259-63/SS 2883 binary sys- 
tem in the IC scattering of soft radiation from the massive 
star. Models with domination of the adiabatic and radiative 
energy losses were considered. However, they predicted the 
maximum of the TeV emission around the periastron pas- 
sage which seems to be inconsistent with the observations. 
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However, very recently Khangulyan et al. (2007) argued that 
under special conditions such a general leptonic model can 
explain the TeV 7-ray light curve observed from this sys- 
tem. Also the hadronic model for production of high energy 
7-rays has been considered by Kawachi et al. (2004) who 
postulate acceleration of hadrons at the termination shock 
in the pulsar wind. In this model the hadrons can pass to 
the dense equatorial wind of the massive star and produce 
pions decaying in 7-rays. The 7-ray light curve expected in 
this model shows two maxima corresponding to the situation 
when the pulsar crosses of the equatorial wind of the massive 
star and the local minimum at the periastron passage. Such 
a 7-ray light curve describes much better the HESS observa- 
tions. The hadronic model has been also recently discussed 
by Neronov & Chernyakova (2007). 

The above mentioned leptonic models have not taken 
into account the influence of the effects connected with an 
aspherical wind of the Be star. Such stars are characterized 
by dense but slow equatorial winds and low density but fast 
poloidal winds. The pulsar orbiting around the massive star 
has to pass trough these different winds. When it is im- 
mersed in the equatorial wind, the shock appears relatively 
close to the neutron star surface. In contrast, in the polar 
wind, the shock extends to relatively close to the massive 
star surface. In this last case, particles accelerated at the 
shock propagate through a much stronger radiation field 
and the ICS process should occur more efficiently. More- 
over, the 7-rays produced close to the stellar surface can 
also be partially absorbed, initiating in this way the IC e 
pair cascade. The importance of the 7-ray production in the 
cascades occurring close to the surface of the massive star 
has been considered in the past, but only in the case of the 
isotropic winds (e.g. Sierpowska & Bednarek 2005a). In this 
paper, we consider those leptonic processes in a much more 
complicated geometry, i.e. by applying the aspherical wind 
properties. We calculate the 7-ray spectra and light curves 
expected in such a more complicated and realistic model. 
The results of our calculations are compared with the recent 
observations of this binary system by the HESS telescopes 
in the TeV 7-ray energies. The preliminary results has been 
published in Sierpowska & Bednarek (2005b) and more de- 
tails of the calculations are included in the PhD thesis of 
Sierpowska-Bartosik (Sierpowska-Bartosik 2006) . 
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Figure 1. Definition of the basic parameters which describe the 
complex geometry of the binary system. iVp denotes direction 
towards the observer, iV is the normal to the orbital plane and i 
is the inclination angle of the orbital plane. The inclination of the 
massive star equatorial wind plane to the orbital plane is denoted 
by while is the normal to the plane of the massive star 
equatorial wind. The longitude of the periastron passage and the 
longitude of the plane of the equatorial wind with respect to the 
periastron passage (defined in the same way) are marked by w p 
and u d , respectively. 
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2 THE BINARY SYSTEM PSR 
B1259-63/SS2883 



The pulsar PSR B1259-63 was discovered during a ra- 
dio survey of southern Galaxy plain in 1992 (Johston et al. 
1992a). It has a period of P = 47ms and period deriva- 
tive P = 2.28 x 10~ 15 . The surface magnetic field has been 
estimated to -Bns ~ 3.3 x 10 11 G based on the above val- 
ues and the rotating magnetic dipole model. The pulsar is 
in an orbit around a massive star of Be type with eccen- 
tricity e = 0.87 and period 1236.72 days (Johnston et al. 
1992b, 1994). The system is viewed by the observer at an 
inclination angle ~ 36° (Wex et al. 1998). The separation 
between the stars changes within the range (0.65 — 9.3) AU 
((0.97 — 14.0) x 10 13 cm). At the periastron the pulsar is at 
the distance of ~23 R*, while at the apastron at ~177 R* 



Table 1. Parameters of the PSR B1259-63/SS2883 binary sys- 
tem. 



(the radius of the Be star is i?* = 4.2 x 10 11 cm). The effec- 
tive temperature of the massive star (B2e) is ~ 2.7 x 10 4 K 
(Tavani et al. 1997). The distance to the binary system is 
not precisely known, as from the photometric method it is 
about d ~1.5 kpc (Johnston et al. 1996), while from the ra- 
dio dispersion measure ~4.5 kpc (Taylor & Cordes 1993). 
The basic parameters of this binary system are listed in Ta- 
ble 1. Due to the presence of an aspherical massive star wind, 
it is expected that the pulsar passes through the equatorial 
wind at phases around ~ (r — 13 ± 5) days before perias- 
tron (noted as r) and after periastron, within ~ (r + 17±5) 
days (Connors 2002). The geometry of the binary system 
and definitions of the basic orbital parameters are shown in 
Fig. 1. 
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2.1 High energy observations 

The binary system PSR B1259-63 was observed in X-rays 
by the ROSAT and ASCA experiments. The measured lu- 
minosity after the periastron passage was L x — (0.8 — 35) x 
10 33 ergs" 1 (assuming d = 2kpc, Cominsky et al. 1994). De- 
tection of X-rays together with the radio emission implied 
that the accretion of matter on the compact object does 
not occur. Observations at the same periastron passage have 
been also performed by the Compton Gamma Ray Observa- 
tory (Grove et al. 1995, Tavani et al. 1996). The OSSE de- 
tector discovered hard X-ray emission with a flat spectrum 
dN^/dE = (2.8 ±0.7) x 10" 3 (£;/100keV)- a cm- 2 s- 1 MeV- 1 
and spectral index a = 1.8 ± 0.6 (Grove et al. 1995). At 
higher energies, only the upper limits could be set (Tavani 
1996). During the last periastron passage the binary system 
has been also observed by the INTEGRAL satellite in en- 
ergy range (3keV - lOMeV) (Shaw et al. 2004). 

The TeV observations by the CANGAROO telescopes 
resulted only with upper limits at the level of 13% of Crab 
flux at ~47 and ~157 days after periastron (Kawachi et al. 
2004). Close to the periastron, the TeV 7-ray emission was 
detected by the HESS telescopes (Aharonian et al. 2004). 
The differential energy spectrum has been described by a 
simple power law, dN/dE — Fo(E /lTeV)~ a , with the spec- 
tral index a = 2.7 ± 0.2 stat and F = (1.3 ± 0.1 stot ) x 
10 _12 cm _2 s _1 . The spectral index was stable with the phase 
of the binary system but the level of emission was clearly 
variable in the range F = (1 — 2) x 10 _12 cm _2 s _1 . The in- 
tegral average photon flux is F(> 380GeV) = (4.0 ± 0.4) x 
10 _12 cm _2 s _1 corresponding to 4.9% of the Crab units. 
The maximum TeV emission was detected after the peri- 
astron at the phase ~ r + 15 days. After that phase, the 
flux decreases up to the phase r + 75 days. HESS observa- 
tions did not cover phases just around periastron passage 
because of bad atmospheric conditions. The average pho- 
ton flux above 380 GeV corresponds to a source luminosity 
L w 8 x 10 32 ergs _1 (at the distance d — 1.5 kpc) which is 
only ~0.1% of the rotational energy loss rate of the pulsar 
PSR B1259-63. 



2.2 The structure of the massive star wind 

The observations of the massive star SS2883 (IR and UV 
energy range) indicate the complex nature of the stellar 
wind as expected in the case of the Be type star winds (Wa- 
ters, 1988). They are composed with an equatorial and polar 
winds characterized by different density and wind velocity. 
The interpretation of the presence of the equatorial wind in 
this binary system is also supported by the radio observa- 
tions. Pulsed radio emission is eclipsed for about 5 weeks 
around periastron (~ r — 16, r + 15 days) and it is depolar- 
ized for about 200 days after the radio eclipse period. These 
observations are interpreted in terms of the model in which 
the plane of the dense equatorial wind is inclined to the or- 
bital plane of the system (Melatos 1995, Wex 1998). Based 
on such a model the pulsar passes twice through the equa- 
torial wind, before the periastron at the phases (t — 13 ± 5) 
days and after it at the phases (r + 7 ± 5) days (Connors 
2002). 

The equatorial stellar wind is dense and has relatively 



low velocity. On the other hand, the polar wind has much 
higher velocity but its density is much lower is respect to 
the equatorial region. The polar wind velocity can be ap- 
proximated by the profile (Casinelli 1979), 

V p (r) = V + (V oo -V )(l-R t ,/r) p , (1) 

where Voc is the terminal velocity of the wind at infinity, 
Vo — 0.01 Voo is the initial velocity, and the index is usually 
taken equal to (3 = 1.5. For the Be type star one may assume 
that V x ~ (1 — 2) x 10 3 kms -1 . The equatorial wind velocity 
is described by, 

V d (r) = Voir/R*) 1 - 25 , (2) 

where the initial velocity is of the order of Vo = lOkms -1 
(see e.g. Lamers 1987). As a result the mass loss rate 
in the polar region, M p , is much lower than that one in 
the equatorial region Ma- For the Be type star, the rela- 
tion between the mass loss rates of these two regions is 
M p = (10~ 2 - 10~ 4 )M d - The interaction of the pulsar wind 
with such a complex stellar wind is additionally complicated 
by the elongated and inclined orbit of the pulsar. 

2.3 Interaction of the pulsar and stellar winds 

The observations of PSR B1259-63/SS2883 indicate that the 
rotational axis of the massive star is inclined to the orbital 
plane of the binary system. Therefore, the pulsar passes 
trough wind regions with different properties. Assuming that 
the pulsar power is constant with the orbital period, it is ex- 
pected that the termination shock is closer to the massive 
star surface when the pulsar is immersed inside the polar 
wind and it is much closer to the pulsar surface when the 
pulsar is inside the equatorial wind. The geometry of the 
termination shock has to be quite complex when the pulsar 
is close to the transition region between these two types of 
winds. Moreover, the wind cavity around the pulsar is not 
axially symmetric with respect to the axis defined by the 
stars in contrary to the case spherically symmetric winds 
(Girard & Wilson 1987). In the case of the spherical winds, 
the geometry of the shock in the binary system with a pulsar 
of power, Lpsr, is determined by the parameter 

V = L PS r/(cMbcVbc), (3) 

where Mb c and Vb c are the massive star wind parameters 
equal to M p or M d and V p and Vd, depending on the type 
of the wind (polar or equatorial, respectively). 

The case of the wind considered here is much more com- 
plicated even if we assume, as the first approximation, that 
the pulsar wind is spherically symmetric while the massive 
star wind has the complex structure discussed above. The 
equatorial wind region can be defined by the opening an- 
gle 9 d , its inclination i d to the line of sight and the longi- 
tude prior to the periastron uj d (with respect to the orbital 
plane). We find the localization of the shock front by com- 
paring the wind pressures given by the velocity and the mass 
loss rates. The geometry and localization of the termination 
shock depend on the pulsar phase and on the localization of 
the stellar equatorial wind with respect to the orbital plane. 
The pulsar phase is denoted as tp (starting from the perias- 
tron) and is given by the orbital parameters, the longitude 
of periastron tu p and the binary inclination i. It is useful to 
introduce the parameter Alu = io d — lu p (see Fig. 1). If the 
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Figure 2. Formation of the termination shock inside the binary with the separation r p = 23-R* for different inclinations of the massive 
star equatorial wind i^. The XZ plane is perpendicular to the massive star equatorial wind plane with the star in the center. The opening 
angle of the equatorial wind is fixed to 9^ = 30° (shaded area). In Fig. (a) the ratio of the mass loss rate of the polar to equatorial 
wind is M p /Md = 10~ 2 , in (b) 10~ 3 , and in (c) 10~ 4 , while the pulsar spin-down luminosity is from Table 1. The termination shock 
formed for different inclination is plotted with specific lines. The corresponding pulsar position is marked by the number: (1, solid line) 
id = 0°, (2, dashed line) i d = 20°, (3, dot-dashed line) i d = 40°, (4, dotted line) i d = 60°, and (5, dot-dot-dot-dashed line) i d = 90°. (d) 
The dependence of the shock structure on the binary separation. The position of the termination shock for the pulsar distance from the 
massive star r = 5 R+ is plotted with the dotted line, for r = 10 R* with dashed line, for r = 20 R* with dot-dashed line and for r = 40 
R* with solid line. Other wind parameters arc: the equatorial wind opening angle 8d = 30°, its inclination id = 20° and the mass loss 
rate M p = lQ- 3 M d - 



pulsar passes through the massive star equatorial plane at 
phases tp = ±90° this parameter is Aco = 0. 

The specific geometry of the termination shock for dif- 
ferent parameters of the stellar wind, but for the same or- 
bital parameters, are shown in Fig. 2. Note, that in the 
region of dense and slow equatorial wind, r\ is lower than 
unity and the shock front is closer to the pulsar (see e.g. the 
case marked by T' in Fig. 2(a, b, c). On the other hand, 
for the polar wind, r\ is larger than unity and the shock 
front is closer to the massive star surface. For that reason, 
when the pulsar is immersed inside the poloidal wind, the 
shock surface can appear quite close to the massive star. 
Leptons accelerated at these parts of the shock propagate in 
the strong radiation field of the massive star and they can 
comptonize stellar photons very efficiently. 



3 MODEL FOR HIGH ENERGY RADIATION 

The complicated shock structures, which can appear in the 
model taking into account the different proprieties of the 
poloidal and equatorial stellar winds, provides a good op- 
portunity for efficient production of high energy radiation. 
However, the presence of the strong (or weak) radiation field 
can influence the conditions for acceleration of leptons at 
the relativistic shock due to the radiation energy losses. All 
these processes, the acceleration, production of photons and 
their possible absorption, should essentially depend on the 
location of the shock, i.e. on the parameters and geometry 
of the binary system. Below we consider these important 
processes. 



3.1 Acceleration of leptons 

Leptons can be already accelerated in the inner pulsar mag- 
netosphere in the electric field induced by rotating mag- 
netic dipole. These leptons initiate cascades in the radiation 
and/or magnetic field as considered in different models (po- 
lar gap e.g. Ruderman & Sutherland 1975, Daugherty & 
Harding 1982, outer gap e.g. Cheng, Ho, Ruderman 1986, 
slot gap, e.g. Arons & Scharlemann 1979). The charged 
products of these cascades escape through the light cylin- 
der. They can be farther accelerated in the pulsar wind zone 
(Kennel & Coroniti 1984). Here we consider two possibilities 
for the spectrum of primary leptons in the pulsar wind zone 
(PWZ) which can represent the limiting cases. In the first 
model (Model A), we assume that leptons inside the PWZ 
have a power-law spectrum with index -1.2 in the energy 
range 500 MeV — 500 GeV as they are injected through the 
light cylinder (see Hibshman & Arons 2001). In the second 
model (model B), leptons are monoenergetic with energies 
E e = 10 7 MeV (Kennel & Coroniti 1984). 

The leptons which reach the shock region can be also 
accelerated in the classical shock acceleration scenario. How- 
ever, the details of this acceleration process are complicated 
since, as we showed above, the shock structure can be quite 
complex. In order to estimate the maximum energies of ac- 
celerated leptons, we have to determine the magnetic field 
strength at different parts of the shock. Note, that the mag- 
netic field strength at different parts of the shock also change 
significantly with the orbital phase. The acceleration time of 
leptons at the specific location at the shock depends on the 
magnetic field strength from the pulsar side. It can be esti- 
mated by, 

race- (Rl/c)(c/V s ) 2 , (4) 
where the Larmor radius of leptons is Rl = emc 2 ^ / B(R a ) « 
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1007r s /a 1/2 cm, B(R B ) = ^B R 3 PS RRlc r ~ 1 ( the ma S" 
netic field in the termination shock is compressed and finally 
B s =3Bp S ii(rs SH ) (Khangoulian & Bogovalov, 2003)), V s 
is the velocity of plasma in the termination shock, and r s is 
the distance from the pulsar to the termination shock. For 
ultra-relativistic shock waves, the velocity of plasma in the 
frame of the shock approaches V s ~ c/3, (Achterberg et al. 
2001). Then, r acc « 9i? L /c « 3 x 10~ 8 7r s /cr 1/2 s. On the 
other hand, leptons lose energy on radiation processes. The 
synchrotron cooling timescale is, 

"T syn = (mc 2 7)/(o- T c£/ B 7 2 ), (5) 

where <tt is the Thomson cross section, Ub = B 2 /(8tv) « 
2.5 x 10 10 B 2 eVcm- 3 = 5.8 x 10 12 cr/r 2 eVcm~ 3 is the en- 
ergy density of the magnetic field, where B is in Gauss. For 
the shock in PSR B1259-63/SS2883, we get approximation 
T S yn ~ 3.3 x 10 6 r 2 /a7 s. The inverse Compton cooling time 
scale in the stellar radiation (in the Thomson regime) can 
be estimated by, 

tic = (mc 2 y)/(a T cU Tad 'Y 2 ), (6) 

where I7 ra d = 4xtsbT+ t b/b b / c ls ^ ne ener gy density of the 
radiation field given by the surface temperature of the mas- 
sive star T*. For a temperature of the massive star equal to 
T* = 2.7 x 10 4 K and at the distance from its surface R B , 
t/rad « 4.7 x 10- 3 T*r- 2 Bc eVcm- 3 , where r s/Be = R s /Rbc 
is the distance from the surface of the massive star to 
the shock. In the Klein-Nishina (KN) regime (7t/kn ^ 
mc 2 /3&bT*) the cooling time can be approximated by ap- 
plying the Lorentz factor of lepton 7t/kn at the transition 
between the T and KN regimes. In the KN regime energy 
losses on ICS depend only logarithmically on the electron 
Lorentz factor, t™(j) w 10 _1 °7 2 Tj C . For the considered 
massive star, we obtain 

^(7) « 7.7 X 10" Vs%e s. (7) 

The maximum energies of leptons accelerated at the shock 
are obtained from comparison of the energy gained by them 
with energy losses. Provided that the synchrotron losses 
dominate, the maximum energies of leptons are, 

7s -r « lovyv- 1 / 4 . (8) 

If the acceleration is balanced by IC losses in the Thom- 
son regime, 7 <C mc 2 /3£;BT*, the maximum Lorentz factors 
of leptons are 

7l c| x « 5.1 x 10 5 iVB.ff 1/4 r B - 1/2 - (9) 

The acceleration process of leptons can be also limited by 
the time they spend in the acceleration region. This time 
scale is defined by the advection process of leptons along the 
termination shock. The advection time can be estimated by 
dividing the characteristic distance scale of the system by 
the velocity of the plasma after the shock (equal to ~ c/3), 

Tadv « R s /V s « 42 r s s. (10) 

However, this last process becomes important only for lep- 
tons with Lorentz factors 7Sv X w 1.4 x lO 9 ^^ (from com- 
parison of Eq. 4 and 10). 

It is usually argued, that, as a result of acceleration at 
the relativistic shock, leptons obtain a power-law spectrum 
with spectral index, a*. For strong relativistic shocks, the 
spectral index is close to, a; ~ 2.2 — 2.3 (see recent models 



by Achterberg et al. (2001) or Ostrowski & Bednarz (2002)). 
The normalization factor in the spectrum depends on the 
acceleration efficiency. In the case of the Crab Nebula it is 
the order of ~ 10%. 

In the regions of the shock which are relatively close to 
the surface of the massive star, the radiation field is strong 
but the magnetic field is relatively weak due to larger dis- 
tance from the neutron star surface. We estimate the maxi- 
mum energies of locally accelerated leptons at the shock by 
taking into account both, their synchrotron and IC energy 
losses as determined above. 



3.2 Escape of gamma-rays 

For some phases of the pulsar in its orbit, leptons can be ac- 
celerated within a few radii from the surface of the massive 
star. TeV 7-rays produced at such distances can be absorbed 
in the radiation field of the massive star. As a result, an 
IC pair cascade in the anisotropic radiation of the mas- 
sive star can develop as considered e.g. by Bednarek (1997) 
or Sierpowska & Bednarek (2005a). The conditions for the 
propagation of leptons and 7-ray photons in the radiation 
field of the massive star can be analyzed from the optical 
depths for leptons on the IC scattering and for 7-rays on e ± 
pair production. To simplify the problem we assume that 
these particles propagate linearly from the place of injection 
at the distance from the center of the massive star, Xi, and 
at fixed angle, 9 obs , measured with respect to the direction 
defined by the centers of the stars. The results of these cal- 
culations are shown in Fig. 3. It is clear, that leptons and 
7-rays (with energies above a few 10 5 MeV) injected within 
the distance of the periastron passage of PSR B1259-63 and 
propagating within 60° towards the massive star should in- 
teract frequently (Fig. 3a shows the results for r = 23i?*). 
The optical depths to ICS increase significantly for smaller 
distances from the massive star (see Fig. 3b). Therefore, 
we conclude that the 7-ray absorption process has to be 
also taken into account even in the case of not so compact 
binary system as PSR B1259-63/SS2883. 

3.3 Gamma-ray production 

Injected leptons start to produce 7-rays (by IC scattering of 
stellar photons) already inside the pulsar wind zone (PWZ), 
after that, during their propagation inside the shock region, 
and finally in the cascade process during propagation of 7- 
rays in the massive star wind region (MSWR). However, the 
relative contribution of these regions to the 7-ray spectrum 
escaping to the observer depends on the distance of the pul- 
sar from the massive star (the phase of the binary system) 
and on the location of the shock. 

Let us consider at first the pulsar wind zone (PWZ). It 
is assumed that leptons propagate rectilinearly in this re- 
gion. Produced by them 7-rays pass through the shock and 
some of them might be absorbed in the region dominated 
by the wind of the massive star. In this region, only lep- 
tons propagating directly towards the observer can produce 
observable 7-rays. 

Leptons which were able to reach the shock region 
are farther accelerated in the shock acceleration mecha- 
nism. They obtain the power law spectra with maximum 
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Figure 3. The optical depths, r, for leptons to Inverse Comp- 
ton Scattering (ICS) and for 7-rays to the pair production 
process in interactions with soft photons from the massive star. 
The surface temperature of the star SS 2883 is = 2.7 X 10 4 K 
and its radius is ij* = 4.2 X 10 11 cm. (a) The optical depths as 
a function of energy of the interacting particle, assuming their 
rectilinear propagation from the injection place, at r ss 23i?*, 
up to infinity (except for the angle 9 = 180° which is the di- 
rection towards the center of the massive star center). Specific 
curves show the results corresponding to different angles of the 
injection: 6 = 60° (dot-dashed line), 9 = 90° (dashed), 9 = 120° 
(dot-dot-dot-dashed), 9 = 150° (solid), 9 = 180° (dotted), (b) 
The optical depths for leptons to ICS process as a function of 
the distance from the surface of the massive star. The energy of 
leptons is fixed to E e = 4 X 10 5 MeV. The angles of injection are 
the following: = 0° (dot-dot-dot-dashed line), 9 = 30° (dotted), 
9 = 60° (dot-dashed), 9 = 90° (dashed), 9 = 120° (thin solid), 
9 = 150° (thick solid). The gray square marked in the both plots 
shows the value of the optical depth equal to unity for leptons 
with energy E e = 4 X 10 5 MeV, injection angle 9 ~ 150°, and at 
a distance r ~ 23-R*. 



energies determined by the escape or radiation processes 
(see Sect. 3.1). We assume that leptons at the shock are 
isotropised in the rest frame of the plasma moving with the 
velocity v&dv — c/3 along the shock surface. These leptons 
produce 7-rays in the IC scattering of stellar radiation. The 
details of this process are quite complicated since the radi- 
ation field is not isotropic with respect to the location of 
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Figure 4. The scenario for the gamma-ray production process 
inside the binary system with complex structure of the termina- 
tion shock. e+e — pairs escaping from the pulsar light cylinder 
propagate through the pulsar wind zone. They are additionally 
accelerated in the pulsar wind termination shock formed by the 
interaction of the pulsar and stellar winds. Gamma-rays are pro- 
duced by leptons due to the Inverse Compton scattering of stellar 
radiation from the massive star. Those photons which propagate 
close to the stellar surface can be absorbed initiating secondary 
generation of pairs. As a result, IC pair cascade can develop 
provided that a part of the shock structure is relatively close to 
the surface of the massive star. 



leptons at the shock. Therefore, the efficiency of 7-ray pro- 
duction and their spectra depend on their escape direction 
(measured with respect to the massive star). Part of pro- 
duced 7-rays (mainly those moving in the general direction 
towards the massive star, i.e. through the massive star wind 
region) can be absorbed initiating the IC pair cascade. 
Also the synchrotron energy losses of these secondary (and 
primary leptons) are taken into account in our calculations. 
The schematic presentation of processes discussed here are 
shown in Fig. 4. As a final result, we obtain the 7-ray spec- 
tra escaping from the binary system at an arbitrary location 
of the observer and different phases of the compact object. 



4 7-RAYS FROM PSR 1259-63/SS 2883 

Low energy observations of the system indicate that the 
plane of the equatorial wind is inclined with respect to the 
orbital plane of the pulsar. This angle is estimated in the 
range id ~ (10°, 40°). The opening angle of the equato- 
rial wind is limited to the range 9j ~ (5°, 40°). Because 
the inclination of the binary system with respect to the ob- 
server is i — 36°, the range of observation angles, which 
determine the IC e pair cascade process, is limited to 
54° < 9 bs < 126°, where 9 bs is the angle between the 
axis defined by the stars and the direction to the observer. 
The parameters of the wind of the massive star SS2883 are 
not well known. Therefore, we have to cover specific range 
of wind velocities and mass loss rates which are expected to 
be typical for the Be type stars. 

The periods when the pulsar passes through the equato- 
rial wind are determined from the radio observations (Con- 
nors et al. 2002) . We analyze different sets of parameters de- 
scribing the geometry of the relative location of the equato- 
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Figure 5. The geometry of the binary system for specific massive star equatorial wind in model I (Au> d = 5°, 6 d = 20°, id = 40°, left 
figure) and model II (Aw d = 20°, 8 d = 20°, id = 30°, right figure). The shaded areas correspond to the phase ranges when the pulsar 
crosses the massive star equatorial wind. 



rial disk with respect to the orbital plane which give similar 
periods of wind disc crossing by the pulsar. The parameters 
which have to be fixed for every model, 9d, id, and co d , arc 
defined in Fig. 1. Observations of the binary system PSR 
B1259-63/SS288 suggest that the plane of the massive star 
equatorial wind is twisted with respect to the orbital plane, 
but the difference Alo = Lo d — to p is rather small, Au ~ 5°. If 
one assumes that the disc is thin, its inclination to the orbital 
plane has to be grater than id > 40°. For our simulations, 
we choose the intermediate opening angle d = 20° but dif- 
ferent disc inclinations id and uJd are investigated. Two sets 
of those parameters, which are generally consistent with the 
observations, have been investigated (see Fig. 5 for details 
of the orbits) : 

• model I: co d = 142°, Alo = 5°, d = 20°, i d = 40°, 

• model II: co d = 158°, Acu = 20°, d = 20°, i d = 30°. 

In the first model, the longitude of the equatorial plane 
with respect to the periastron phase is fixed to Au) = 5°. 
To achieve the agreement with periods of pulsar passage 
through the equatorial wind, the inclination angle has to be 
equal to i d = 40° (see Table 2). The second model represents 
a more tilted disc, Alo = 20°. In this case, in order to cover 
the observed periods of the disk crossings, the inclination 
angle has to be i d = 30°. 



4.1 Gamma-rays from the pulsar wind region 

As we already noted, leptons, injected through the pulsar 
light cylinder, move radially in the pulsar wind zone. They 
remain almost at rest with respect to the pulsar wind. There- 
fore, their synchrotron energy losses can be neglected at first 
approximation. However, these leptons move in the radiation 
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Table 2. The relation between the orbital phase and the time 
(in days) from the periastron (denoted by r). t ln and t ou t arc 
the time periods corresponding to the pulsar passage through 
the equatorial wind region (Connors et al. 2002) for considered 
models I and II. 



field of the massive star. Thus, they can produce 7-rays in 
the ICS process. We calculate these 7-ray spectra for two 
mentioned above models for the spectra of leptons and ge- 
ometries of the shock (see Sect. 3.1). Note that 7-rays pro- 
duced by leptons in the PWZ can be absorbed in this same 
soft radiation field during their propagation in the pulsar 
wind zone and the massive star wind region. These cascad- 
ing effects are taken into account in our model as described 
in detail in Sierpowska & Bednarek (2005a). The specific 7- 
ray spectra escaping from the binary system were calculated 
for the pulsar at the periastron passage, r = 23ii*. At this 
location, the observer views the binary system at the an- 
gle 0*fa w 113°. The distance from the pulsar to the shock 
front was calculated in each model: n = 13 7?* (model I), 
r2 = 8.8 (model II). As expected from the analysis of 
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Figure 6. The differential 7-ray spectra (SED) produced in the 
pulsar wind zone (PWZ) for the pulsar at the periastron passage 
(r = 23-R*). The angle to the observer is 9 ob3 = 113°. The spectra 
are calculated for leptons with the mono-energetic spectra (energy 
E e = 10 7 MeV, thin lines) and with the power law spectrum 
(see the text, thick lines). The distance from the pulsar to the 
termination shock is ri = 13 R* (solid line), T2 = 8.8 R* (dashed 
line) for the geometrical models I and II. 



the optical depth, the 7-ray fluxes are larger for larger vol- 
ume of the PWZ (the shock is farther from the pulsar). For 
monoenergetic primary leptons (model A), the spectrum of 
7-rays is very flat (spectral index ~ 0.5). In case of injec- 
tion of primary leptons with the power law spectra (model 
B), 7-ray spectrum shows maximum corresponding to IC 
scattering of stellar photons in the KN regime. The gen- 
eral spectrum is significantly steeper in this case (spectral 
index ~ 1.5 below ~ 3 x 10 4 MeV). Similar spectra have 
been obtained in the case of monoenergetic leptons inside 
the pulsar wind nebulae (see Ball & Kirk 2000, Bogovalov 
& Aharonian 2000, Khangulyan et al. 2007). Observation of 
such very flat spectra from the binary system would indicate 
the efficient acceleration of leptons already inside the pulsar 
magnetosphere or in the vicinity of the pulsar light cylinder. 
In the case of PSR B1259/SS2883 such flat gamma-ray spec- 
tra have not been detected up to now.Therefore, the PWZ 
seems not to be responsible for the bulk of 7-rays escaping 
from this binary system. However, in other binaries, in which 
the pulsar wind region completely dominates the volume of 
the binary system, e.g. as in the case of the famous mil- 
lisecond binary system PSR 1957+20, such very flat 7-ray 
spectra might be produced. Their observation should give 
important information on the spectra of primary leptons in 
the pulsar wind zone. 

4.2 Gamma-rays from the shock region 

Since the radiation processes of leptons, accelerated at the 
shock region, depend significantly on the details of the shock 
geometry, we calculate the 7-ray spectra for two models 
which fulfill the general constraints concerning the positions 
of the pulsar close to the periastron passage. We fixed the 
basic parameters characterizing the wind of the massive star: 
the equatorial wind mass loss rate is M d = 2.0 x 1O~ 7 M 
year -1 , the polar wind M p = 10~ 3 Md, i.e consistent with 
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Figure 7. The geometry of the termination shock with respect 
to the massive star equatorial wind (the massive star in the cen- 
ter) for the subsequent pulsar phases in the model I and II. The 
pulsar phases are denoted with numbers while the shock shapes 
are plotted with different line styles: (1) <p = —90° - thin dotted 
line, (2) <p = -60° - dotted one, (3) tp = -30° - dot-dot-dot- 
dashed, (4) ip = 0° - solid, (5) <p = 30° - dashed, (6) ip = 60° - 
dot-dashed, (7) ip = 90° - thin dot-dashed. 



the observations of this system. Other parameters of the sys- 
tem derived from different observations are shown in Table 
1 . As the orbit is highly eccentric and the plane of the equa- 
torial wind is inclined to the orbital plane, the shape of the 
shock structure changes from phase to phase in a compli- 
cated way. 

We apply two geometrical models mentioned above. 
The details of these models are shown in Fig. 5 and the 
shapes of the shock front in subsequent phases of the bi- 
nary system are shown in Fig. 7. Note that for three shown 
phases, ip = —90°, —60°, 90°, the pulsar is immersed in- 
side the equatorial wind. For both models, we calculate the 
7-ray spectra escaping to the observer by applying the modi- 
fied Monte Carlo code developed by Sierpowska & Bednarek 
(2005a). The spectra of primary leptons injected at the shock 
are normalized to one particle per sr. Calculated 7-ray spec- 
tra allows us to obtain the 7-ray light curves in two energy 
ranges which correspond to the energy ranges of the GLAST 
detector (1 GeV— 100 GeV) and the present Cherenkov tele- 
scopes (> 100 GeV). We also investigate the 7-ray spectra 
and light curves for different power law spectra of primary 
leptons applying the range of spectral indexes oti = 2.0 — 3.5. 
Calculated phase dependent 7-ray spectra are shown for the 
spectral index of primary leptons a; = 2.5 (Fig. 8). As ex- 
pected, the spectral index of produced 7-rays is equal to 
q 7 ~ (2.2, 2.3) for energies < lOOGeV. It is clearly related 
to the spectral index of primary leptons when the scattering 
occurs in the Thomson regime. 7-ray photons with larger en- 
ergies are created in the Klein-Nishina regime, and so their 
spectra steepens. Note, that for large range of phases around 
periastron, the spectrum of 7-rays, calculated in terms of the 
model I, does not change drastically. However, if the pulsar 
is farther from the periastron, the cut-off in the 7-ray spec- 
trum appears at lower energies. Note that calculations of the 
7-rays spectra in such a complicated scenarios (which addi- 
tionally take into account the cascade processes) are very 
time consuming. Therefore, the spectra suffer to some ex- 
tent from low statistics. However, their main features can 
be clearly identified. 
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Figure 8. The differential 7-ray spectra (SED) for subsequent 
pulsar phases in the case of model I parameters (figures (a) and 
(b)) and model II ((c) and (d)). Primary leptons are injected 
with the power-law spectra and spectral index oti = 2.5. The 7- 
ray spectra are shown for different pulsar phases before pcriastron 
in figures (a) and (c): tp = —90° (dotted line), ip = —60° (dot- 
dashed), ip = —30° (dashed), tp = 0° (solid gray line), and after 
periastron in figures (b) and (d): (p = 30° (dashed line), ip = 
60° (dot-dashed), <p = 90° (dotted) and ip = 105° (dot-dot-dot- 
dashed). The magnetization parameter of the pulsar wind is fixed 
to a = 0.01. 



In Fig. 9 we also show the 7-ray light curves in the 
GLAST and Cherenkov telescopes energy ranges. The 7-ray 
light curves calculated in terms of the model I are charac- 
terized by clear maximum at the periastron passage and the 
second maximum at the phase r + (15— 17) days, (<p ~ 90°), 
corresponding to the second passage of the pulsar through 
the equatorial wind of the massive star. The 7-ray light 
curves obtained in terms of the model II are different. They 
do not show the maximum at the periastron passage but 
only slight increase of emission at the second passage of the 
pulsar through the equatorial wind. 

Note that the equatorial wind of the massive star is ori- 
ented with respect to the observer in a different way in both 
discussed models. The 7-rays, produced in the shock region, 
propagate along different paths to the observer (i.e. closer 
or further to the massive star), but absorption of photons 
can occur effectively only in the first case. The differences 
can be observed when we compare the photon spectra be- 
fore and after periastron. Moreover, 7-ray fluxes are higher 
after periastron phase which is due to the fact that pulsar 
changes its position with respect to the observer after pe- 
riastron (begins to approach the observer and the massive 
star is 'behind' it). Also the orientation of the shock plane 
changes with respect to the observer. Leptons, which prop- 
agate along the shock, are moving in direction closer to the 
observer line of sight. So then, produced 7-ray spectra can 
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Figure 9. The gamma light curves in two energy ranges: 1 GeV 
- 100 GeV (a, c) and above 100 GeV (b, d). The results for the 
model I are presented in figures (a) and (b) while for the model II 
in (c) and (d). The leptons were injected with the power law spec- 
trum and spectral index ai = 2.5 (the magnetization parameter 
is a = 0.01). The 7-ray luminosity, Lbwi i s in units of photons 
MeV sr™ 1 . The phases corresponding to the pulsar in the massive 
star equatorial wind are shaded. 



be observed without significant absorption having higher en- 
ergies (the optical depths for pair production are the lowest 
when photons propagate outward the surface of the massive 
star). 



5 COMPARISON WITH OBSERVATIONS OF 
PSR B1259-63/SS2883 

In the above discussed models I and II, we applied the 
value for the pulsar wind magnetization parameter equal to 
a = 1CT 2 (as it was done in previous modelling, see Intro- 
duction). However, in such a case the energies of produced 
7-rays are too low in order to be consistent with those ob- 
served from the binary system PSR B1259/SS2883 (Aharo- 
nian et al. 2005). The 7-rays with larger energies can ap- 
pear for the case of the lower magnetization parameter of 
the pulsar wind. Then, due to the weaker magnetic field 
at the shock region, leptons can be accelerated to higher 
maximum energies. Therefore, when comparing the result 
of the model with the observations of PSR 1259-63/SS2883, 
we apply the value a — 10~ 4 . Moreover, we assume that the 
energy spectra of primary leptons have low energy cut-offs 
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Figure 10. The comparison of the 7-ray light curves obtained in 
the model II (a) and model I (b) with the observations of PSR 
B1259-63/SS2883 at energies > 380 GeV (Aharonian et al. 2005). 
The primary spectrum of leptons has the spectral index 2.5 and 
the magnetization parameter of the pulsar wind is a = 10~ 4 . 
The periods when the pulsar crosses the massive star equatorial 
wind are shaded. The 7-ray light curve obtained in model I but 
for different orientation of the orbital plane and the plane of the 
equatorial wind (with uj v = 148°) is shown in (c). 



at E m i n — 10 MeV. We suggest that this low energy cut- 
off corresponds to the energies of leptons which arrive to 
the shock region from the pulsar. The gamma spectra and 
light curves are normalized to the fluxes observed from the 
binary system PSR B1259-63/SS2883 by the HESS experi- 
ment (Figs. 10 and 11). The 7-ray light curve in the energy 
range >380 GeV has two distinct maxima at phases corre- 
sponding to the pulsar crossing the equatorial wind. Max- 
imum 7-ray emission appears at phases symmetrical with 
respect to the periastron, r ± 17 days (p = ±90°). The 7- 
ray light curve for model II (Fig. 10a) has local maximum 
at the periastron passage and the second, stronger one, after 
periastron at phase ~ r + 8 days. Before the periastron the 
data are well fitted in case of both models. For the investi- 
gated model I (Fig. 10b) we get a good general agreement 
with the HESS results, e.g. the power decreasing in period, 
r — 11, r — 5 days, and increasing in, r + 11, r + 20 days. 
A still better fit is obtained when we change by ~ 10° the 
angle lo p which defines the position of the periastron with 
respect to the observer (Fig. 10c). Because the observation 
did not cover the period ~ r — 5, r + 10 wc cannot definitely 
exclude that model. 
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Figure 11. The differential 7-ray spectra (SED) for subsequent 
pulsar phases calculated in model I (with a = 10 -4 ). The 7-ray 
spectrum from the binary system PSR B1259-63/SS2883 observed 
by the HESS telescopes is plotted for comparison (with the errors 
marked by the gray lines). The spectral index of primary leptons 
is oti = 2.5. The spectra are calculated for the phases before pe- 
riastron (a): tp = —90° (dotted line), ip = —60° (dot-dashed), 
<p = —30° (dashed), <p = 0° (solid red line), and after periastron 
(b): tp = 30° (dashed line), p = 60° (dot-dashed), p = 90° (dot- 
ted), ip = 105 (dot-dot-calculated)). 7-ray spectra calculated in 
terms of the model II are shown in (c) and (d). 



The photon spectra for both models and a = 10~ 4 are 
compared to the HESS spectrum in Fig. 11. Photon fluxes 
at phases after periastron (ip > 0) are larger than before 
it. The photon spectra were calculated for different initial 
power law spectra of lepton. The best description of the 7- 
ray spectrum measured by HESS is achieved in the case of 
leptons injected with the power law spectrum and spectral 
index ai — 2.0 — 2.5. Note that for these parameters 7-ray 
spectra can extend up to ~ lOTeV. For index oii = 2.0 and 
energy range 5.5 <log(E 7 )< 6.5 MeV, we get the photon 
index is a 7 ~ 2.5 while for a; = 3.5 the photon spectrum is 
steeper (a 7 ~ 4.0). 



6 DISCUSSION AND CONCLUSION 

We considered the model for production of high energy 
gamma rays in the binary system consisting of a Be type 
massive star and the energetic pulsar applying the details of 
complicated geometry of the massive star wind. The model 
predictions are compared with recent observations of the 
binary system PSR B1259-63/SS2883. In this system, due 
to the interaction of the pulsar wind and the massive star 
wind, a shock wave is created. The geometry of the shock 
is complex and changes with the orbital phase as the pulsar 
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passes through equatorial or polar wind of the massive star. 
Such structure of the massive star wind is characteristic for 
a Be type star. In the considered system, the equatorial wind 
plane is inclined to the orbital one. 7-rays are produced by 
leptons which are accelerated at the shock region and after 
that scatter thermal photons from the massive companion. 
Since the soft photon field is anisotropic, we take into ac- 
count the details of the geometry of the binary system. The 
change of the separation between the stars with the phase 
of the binary and the complex structure of the massive star 
wind have also big impact on the conditions for acceleration 
of leptons and their interactions with the stellar radiation. 
The maximum energies of leptons and their energy losses 
are determined by the magnetic field strength at the shock 
region which in turn is defined by the magnetization pa- 
rameter (a) of the pulsar wind. In our calculations a value 
of a = 1CP 2 and 10~ 4 is considered. The presented model 
takes into account the complicated geometry of the system 
and still it has to depend on many parameters which de- 
scribe the position of the orbital plane and the equatorial 
plane of the stellar wind with respect to the observer. 

In such a complicated geometrically model, leptons can 
be accelerated in some cases relatively close to the stellar 
surface. Therefore, the possible absorption of the 7-rays pro- 
duced in the IC scattering process has to be also taken into 
account. The optical depths for leptons on ICS within the 
pulsar wind region are clearly too low for effective produc- 
tion of 7-ray photons. However, leptons, which arrive from 
the pulsar to the shock region, are trapped by turbulent 
magnetic field after the shock and can be additionally accel- 
erated. If the shock region is relatively close to the massive 
star, IC scattering becomes important and high energy 7- 
rays are produced effectively. 

We calculate the 7-ray spectra produced in the shock 
region for power law spectra of leptons in two different ge- 
ometrical models. By applying the sets of parameters, I: 
Acu d = 5°, 9 d = 20°, i d = 40° and II: Auj d = 20°, 6 d = 
20°, i d — 30°, we reach the general consistency with the TeV 
7-ray observations. The best fit of the gamma light curve 
measured by the HESS telescopes (Aharonian, et al. 2005) 
is obtained by applying the model I with a = 10~ 4 . How- 
ever, this fit has been reached assuming that the periastron 
position of the pulsar is defined by the angle ui p = 148° con- 
trary to the present estimates, i.e. lu p ~ 138°. Note, however 
that the estimation of the angle cu p depends on massive star 
wind parameters such as the inclination of the disc (applied 
from Wex et al. 1998) and the mass loss rate, which in case 
of PSR B1259-63 was fixed to 5 x 10~ 8 M Q yr _1 , and which 
are in fact not very well known. 

The calculated 7-ray spectra are in agreement with the 
TeV HESS data provided that the spectral index of the spec- 
trum of leptons is a* = —2.5. The slope and also the fluxes of 
the photon spectra do not change with orbital phase, as was 
reported by the HESS collaboration. The final light curve 
for model I (a = 10~ 4 , u = 148°) is characterized by a max- 
imum for the phases when the pulsar crosses the equatorial 
disc after periastron. This fit differs from the observational 
data for phases when pulsar approaches the equatorial disc. 
This can be due to the complex geometry and structure of 
the shock region when the pulsar wind interacts with the 
strong equatorial wind. The model assumes that the accel- 
eration of leptons occurs locally at the shock. In reality, this 



acceleration process of leptons arriving to the specific place 
at the shock can extend to other parts of the shock. Such 
a non-local acceleration process is very difficult to consider, 
but it might be important and influence some details of the 
calculated 7-ray light curves. 

Based on the reported calculations, we conclude that 
the general shapes of the TeV and GeV 7-ray light curves 
close to the periastron should be quite similar. However, 
farther out from the periastron, the GeV emission should 
continue with larger flux. This feature can be tested by the 
future observations in the GeV energy range by the AGILE 
and GLAST telescopes. 
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